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In the Terahertz (THz) domain, we investigate both numerically and experimen- 
tally the directional emission of whispering gallery mode resonators that are per- 
turbed by a small scatterer in the vicinity of the resonators rim. We determine 
quality factor degradation, the modal structure and the emission direction for var- 
ious geometries. We find that scatterers do allow for directional emission without 
destroying the resonator's quality factor. This finding allows for new geometries and 
outcoupling scenarios for active whispering gallery mode structures such as quantum 
cascade lasers and passive resonators such as evanescent sensors. The experimental 
results agree well with finite difference time domain simulations. 
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I. INTRODUCTION 

The THz frequency range, denned as the frequency range from 100 GHz to 10 THz, cor- 
responding to a wavelength between 3 mm and 30 fim offers unique advantages for studying 
whispering gallery mode (WGM) resonators. On the one hand, much larger resonators and 
feature sizes have to be used compared to the optical domain. In the lower THz range 
around 100-300 GHz, the wavelength is roughly a factor of 3000 larger than in the visible 
domain. However, optical and physical effects remain the same due to the scalability of 
electromagnetism. It is therefore much easier to characterize the near and far field [8j. On 
the other hand, materials are less advanced than in the optical domain, resulting in losses 
and therefore much lower quality factors. Furthermore, there is a lack of powerful, but at 
the same time, tunable continuous-wave sources. We developed a THz system based on 
n-i-pn-i-p superlattice photomixers [10] that allows for characterizing WGM resonators 
between 60 GHz and several 100 GHz [8]. 

In this paper, we study the emission characteristics of WGM resonators in the lower THz 
range. For most applications of WGM resonators, not only the modal distribution within the 
resonator, but also knowledge on the radiative losses due to the shape and structure of the 
resonator are extremely important. Several detection schemes employ the evanescent field 
surrounding a WGM resonator for sensing. Gas, fluids, or particles surrounding the resonator 
influence the evanescent field by either shifting the effective refractive index or by absorption 
[5j [20]. Consequently, this alters the modal structure of the resonator, resulting in mode 
shifts and altered quality (Q) factors. For WGM lasers, particularly THz quantum cascade 
lasers [I], the directional out-coupled emission is paramount. Therefore, the resonator shape 
often deviates from the circular geometry. Such examples include chaotic micro resonators [2] 
where the light is confined on a stable periodic lightray-orbit or by the short term dynamics 
of the chaotic light ray-dynamics inside the resonator which channels a direct emission port 
[T5| [IB"] . We present another way to achieve directional emission: In this paper, we study 
directional emission from a WGM resonator with a finite scatterer both experimentally and 
theoretically in the THz frequency range. Theoretical calculations by Wiersig et al. [22] 
already described efficient outcoupling of a high Q mode via interaction with a low Q mode 
without ruining the Q factor. Here, however, we address the mode directly by perturbing 
the mode with the scatterer. We will show that the perturbation may be weak enough to 
maintain a considerably high Q factor while the power is coupled out directionally. We will 
use passive resonators that are coupled to a waveguide. The paper is structured as follows: 
First, we measure the influence of the waveguide on the modal structure and near field of 
the resonator by a near field probe. Second, we experimentally determine the influence of 
the position of a hole acting as scatterer on the Q factor of the disk. Third, we will explain 
the numerical algorithm that was used to theoretically determine the radiation pattern of a 
disk with hole. Last, we will show experimental results on the outcoupling performance and 
compare the results to theory. 

II. EXPERIMENTAL SETUP 

We use a n-i-pn-i-p superlattice photomixer [9] as THz source that has been developed in 
house. The measurement setup is illustrated in Fig. [I]a) and b). The source is mounted on a 
silicon lens for pre-collimation and efficient out-coupling of the THz beam. The beam is then 
collimated with a parabolic mirror (PM1). A second parabolic mirror (PM2) focuses the 
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Figure 1. Schematic setup. The resonator is mounted at a fixed position in the center of a rotation 
stage. The detector (Golay cell) scans the angular far held pattern of the disk. Two cylindrical 
lenses (horizontal lens, HD, and detector lens, DL) improve the signal to noise ratio by focusing 
the emitted power from the resonator on the detector. For most measurements, the lenses collect 
radiation within a 10° angle. Alternatively, the lenses can be replaced with a probe waveguide that 
touches the boundary of the disk to scan the near held, b) Photograph of the resonator in the 
center of the ration mount, c) Definition of variables used throughout the paper. Angles are dehned 
with respect to the coupling position, the angle of incidence is dehned as positive in the direction 
of the light. The thickness of the boundary between hole and disk circumference is labelled as b. 
The hole (radius 77J position relative to the coupling position is dehned as ah- Parameters for the 
numerical calculation, such as the incidence angle of a light ray on the surface, x> and the far held 
emission angle 9 are also included. 

THz beam on a feed horn (FH) with a rectangular Teflon waveguide (dimensions 1 mm x 1.5 
mm) mounted inside the horn. The feed horn is used to increase the coupling efficiency. The 
Teflon waveguide delivers the THz power to the WGM resonator. A Golay cell detector (D) 
is used to monitor the transmission through the waveguide. The WGM resonator is situated 
in the center of a rotation stage. The Golay cell can be mounted on the rotation stage to 
scan the far field emission from the resonator. A horizontal cylindrical (HL) lens collects 
THz power emitted perpendicular to the detection plane. A second, vertical cylindrical lens 
(DL) is mounted in front of the Golay cell to collect power within a 10°-15°angle to improve 
the signal to noise ratio. The scan range was limited by the Teflon waveguide and the size 
of the Golay cell and the lenses to about 140°. Alternatively, the Golay cell detector can 
also be attached to a probe waveguide that touches the resonator in order to scan the near 
field. In order to reduce the necessary integration times of the Golay cell novel field-effect 
transistors could also be used pj], [12] 

The frequency resolution (120 MHz) of the system is sufficient to measure Q factors up to 
1200 in the frequency range from 100 — 300 GHz. Due to the lack of ultra-low-loss materials 
at THz frequencies, this resolution is sufficient since the ideal loss-limited Q factors are in 
the range of a few 1000. As an example, the imaginary part of the dielectric constant of 
polyethylene and Teflon [6] is in the range of e" > 1.8 x 10~ 3 , resulting in an absorption- 
limited Q factor of < 1100, even if no radiative losses are present. For the investigated 
disks and frequency ranges, however, the measured Q factor was always smaller than 1200, 
limited by radiative losses. 

In a first step, we investigate the perturbation of the coupler on the modal structure of 
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a) Transmission b) Outcoupled 
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Figure 2. The evanescent near field of a polyethylene WGM resonator (diameter 25 mm, in the 
following referenced as "large resonator" ) without hole is probed with a waveguide, a) Transmission 
through the feed waveguide while the near field probe is scanned along the circumference of the 
resonator, b) Outcoupled power through the near field probe. 



the resonator without scatterer. An unperturbed circular resonator should show a position- 
independent modal field strength. We resolved the local field strength by a near field tapered 
probe waveguide that weakly couples to the disk. In Fig. [2] a) we show the transmitted power 
through the feed waveguide and in Fig. [2] b) the outcoupled power through the probe while 
scanning both frequency across a resonance and the position of the probe. As expected, 
the transmitted power through the feed waveguide is high, when little power is coupled 
out through the probe waveguide. However, we also see a strong, periodic probe position 
dependence of the power in both sub-figures with a five fold symmetry. This number is 
much too small to be the modal number (~ 50) for a resonator with a diameter of 25 mm at 
a wavelength of 2.24 mm. Already the coupling to the feed waveguide perturbs the modal 
structure of the resonator considerably. The orbit of the mode is not circular any more, there 
exist positions where the mode is close to the surface, allowing for enhanced outcoupling, 
and further away, showing local minima in Fig. |2]b). 

As a next step, we investigated the influence of the hole position on the quality factor 
of the resonances with disks with a diameter of d = 10.5 mm at wavelengths between 2.3 
mm and 0.91 mm (0.13-0.33 THz). For holes very far away from the circumference, the 
outermost radial modes that are excited by the waveguide, are not affected by the hole 
as their mode volume does not overlap with it. The Q factors are similar to that of a 
resonator without hole. When holes are drilled closer to the circumference, the overall Q- 
factor decreases as shown in Fig. [3] and the mode position shifts because the mode gets 
perturbed by the hole. The effect becomes more pronounced if the boundary thickness b 
(see Fig. [T] c)) becomes comparable to or smaller than the wavelength. The loss-limited 
quality factor can be estimated using an empiric fit of the form 

Qioss + Qo exp (71/) ' 

where 7 and Qo are fitting factors describing the increase of Q with increasing frequency. 
The formula takes the roughly exponential increase of the quality factor of an unperturbed 



5 




Frequency / GHz Frequency / GHz 

Figure 3. a) Q factors of several disks without and with a hole with a boundary thickness of b 
(see Fig. [I]). The closer the hole is drilled to the resonator circumference, the stronger are the 
perturbations of the mode and the stronger the reduction of the Q factor. The inset shows the 
extrapolated loss-limited Q factor in the limit of large frequencies, b) Modal structure of the disk 
with b = 0.38 mm before and after the hole was drilled. Despite a strong shift in the resonance 
frequency, the Q factor was only slightly altered. 

resonator with increasing frequency into account, that is finally limited by losses to a max- 
imum value of Qi OS s- Qo <C Qioss is the extrapolated quality factor at low frequencies. The 
inset of Fig. [3] shows the extrapolated values for Qi OS s- The disks with b = 0.38 mm and the 
disk with the slit are already loss-limited within the measurement range. 

The Q factors did not degrade drastically due to the perturbation of the hole. This 
demonstrates that finite scatterers can be used for resonators and laser cavities. In the 
following, we discuss the theoretical framework for determining the outcoupling performance 
by the hole. 

III. THEORETICAL MODELING 

We now turn to our theoretical approach to model the experimental data. Such modeling 
is desirable mainly for two reasons. On the one hand, the model gives insight into the 
physical mechanism for the directional emission, not least since it provides access to system 
properties which may be difficult to obtain in the experiment. On the other hand, in 
order to design the emission pattern of WGM resonators towards specific applications, an 
accurate theoretical modeling is desirable in order to perform numerical parameter studies. 
For stationary resonance problems a number of numerical methods exist, such as scattering 
methods [TH], boundary integral methods jT8l ED, E3] and multipole methods [TS] . These 
methods however do not allow for studying the effects of the coupled waveguide. 

The light propagation is thus simulated by the Finite-Difference Time-Domain (FDTD) 
Method (3j Uni E]. The geometry of waveguide and resonator as well as the material 
properties are chosen as in the experiment, without fine-tuning of any of the parameters. We 
however restrict the simulation to a two dimensional rectangular grid in the resonator plane, 
with grid constant Ax = Ay = 0.02 mm. This corresponds to about A/90 or A/110 for the 
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two wavelengths in the mm range we consider below. In two dimensions, Maxwell's equations 
for the six components of E and H separate into two disjunct sets of three equations each. 
As in the experiment, we study the propagation of the EM field in TE 2 polarization, such 
that H z , E x and E y need to be considered. We drive the H z component using a soft source 
placed at one edge of the waveguide, with a transverse mode profile obtained by solving 
the corresponding Helmholtz equation. Unwanted reflections from the boundaries of the 
simulation area are suppressed by Berenger type perfectly matched absorbing boundaries [1] . 

To predict the emission pattern, we proceed in three steps. First, we excite the system 
with a temporally short and spectrally broad pulse and monitor the transmission at the 
waveguide edge opposite to the source. By relating this result to a reference calculation 
without the resonator, we obtain the modal spectrum of the resonator. We then choose a 
single resonator mode by matching the numerically obtained mode spectrum to the experi- 
mentally observed one. 

In the second step, we excite the system with a monochromatic continuous-wave field at 
the resonance frequency of the mode identified in the first step, and evolve the system into its 
stationary state (modulo the oscillations at the incident field frequency). This evolution can 
be monitored, e.g., via the time-dependent power flux through the waveguide downstream 
the resonator. The steady state, where the resonator is fully loaded, is reached on time 
scales short enough to neglect the initial cavity build-up time. 

In the third step, we calculate the emission pattern of the resonator. A principle problem 
arises from restraining the computation time to a moderate level. The FDTD simulations 
therefore are constrained to a small region around the resonator (~ tens of mm) which are 
small compared to the experimental detector distance (~ 10 cm). Thus it is impractical to 
include the detector itself into the simulation. But while the detector exclusively monitors 
the far field radiation component of the emission, the vicinity of the resonator is strongly 
influenced by near field contributions. Thus, a transformation of the FDTD results into the 
far field is required [31 [13J [17] . To this end, we calculate the amplitude and the phase of the 
electric and magnetic fields on a circle around the resonator center, with radius larger than 
the resonator radius by a variable offset p. At each FDTD grid point f on the circle, we fit 
the FDTD time evolution over few cycles of the incident field to the function a sin(wt + </>). 
Here, u is the frequency of the incident field; the amplitude a > and the phase (p are the 
fit parameters. From the field configuration, we evaluate the time averaged Poynting vector 
S(f) at each point f*on the circle. Next, we determine the intersection point R(r) of a ray 
starting at r in the direction S with the detection sphere of radius R around the resonator 
center. As a result, the point r contributes with magnitude \S\ to the detection signal at 
detector position R(r). By incoherently summing up the contribution of all points on the 
circle around the resonator we obtain a prediction for the emission pattern as observed by 
the detector. 

We compared the emission patterns evaluated for different offsets p from the resonator 
edge, and found that for small offsets, the results strongly depend on p. However, already 
with offsets p of few wavelengths of the incident light field, the radiation pattern becomes 
largely independent of the offset, indicating the suppression of near field components. It 
should be noted, however, that neither the simulation range nor the detector distance in the 
experiment are large enough to neglect the finite size of the resonator. 
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Figure 4. Small Resonator (R = 5.2 mm) with a hole at 180° ±5°. a) Snapshot of the FDTD 
calculation on a logarithmic scale. The small circle indicates the hole. b)Result of the Poynting 
vector analysis. The background shows a snapshot of the FDTD simulation. The arrows are 
the Poynting vectors, originating from the point at which they are evaluated, respectively. The 
Poynting vectors clearly reflect the pattern of the radiated field. The distance between the resonator 
edge and the circle on which the Poynting vectors are calculated is p = 8 mm. c) Experimental data 
(black dots) as well as far-field radiation pattern predicted from the theoretical Poynting-analysis 
(red) on a linear scale. The dip in the theoretical emission intensity at 270° is due to the waveguide 
region which was excluded from the Poynting vector analysis. In the experiment, the waveguide is 
curved and thus does not lie in this direction. 



IV. COMPARISON OF THEORETICAL AND EXPERIMENTAL RESULTS 



We investigated two types of resonators for comparison with our theoretical framework. 
First we consider a polyethylene WGM resonator with R = 5.2 mm, = 0.51 mm, b = 
0.41 mm and = 90°. The relative permittivity is e-wgm = 2.56 {n PE = 1.6) for the 
resonator and Eh = 1 (air) inside the hole. This resonator is referred to as small resonator 
in the following. Its resonance frequencies and Q factors are illustrated in Fig. [3] The far 
field emission was characterized for angles between 160° and 290° at a resonance frequency 
of 170 GHz (A = 1.76 mm). 

The dimensions of our second resonator are R = 12.5 mm, = 1.25 mm and b = 1 mm. 
In the following this resonator is referred to as large resonator. The larger size was chosen 
to have less radiative losses of the unperturbed resonator and to simplify probing the near 
field with a probe waveguide. The large resonator was studied at cut = 47° and ah = 315°. 

Fig. [4] a) shows a FDTD snapshot indicating the considered geometry. Fig. [4] b) depicts 
the result of the Poynting vector analysis for the small resonator. The Poynting vector 
structure clearly resembles the interference pattern visible in the field configuration. How- 
ever, the radiation pattern cannot be read off directly from the Poynting vector structure, 
as the projection of the pattern onto the detection sphere has not yet been applied. The 
experimental data are compared to the prediction for the far-field pattern from the Poynting 
analysis in Fig. |4]c). The comparison between experiment and theory in (c) shows excellent 
agreement. The theoretical curve was obtained for a detector size of 10°, and a detector 
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distance of 10 cm, consistent with the experiment. Apart from an overall scaling, no further 
free parameter had to be adjusted. The scaling was chosen such that the largest relative 
emission from the theoretical and the experimental data coincide. The main lobes as well 
as the intensity minima observed in the experiment are clearly reproduced in the theoretical 
analysis. The main difference arises at 270°, where the theoretical data exhibits a reduction 
in the emission intensity. This most likely is due to the waveguide, which in the theoretical 
calculation is situated at this emission angle, such that an angular range around 270° had 
to be excluded from the Poynting vector analysis. In the experiment, the waveguide was 
bent, such that it did not interfere with the emitted radiation at this angle, see Fig [TJ The 
Poynting vector analysis and the comparison between experiment and theory for the large 
resonator with hole at ah = 315° are shown in Fig. |5j Again, a detector size of 10°, and a 
detector distance of 10 cm were chosen. Good qualitative matching between experiment and 
theory is achieved, even though the quantitative agreement is not as good as for the small 
resonator. The main difference is the width of the second lobe around 130° — 180°, which 
is wider in the experimental data than in the theoretical prediction. One possible origin 
for this could be slightly incorrect values for the index of refraction or the distance to the 
coupling waveguide to the resonator which could not be measured accurately, altering the 
perturbation of the resonator. Furthermore, the disk center is centered in the rotation stage. 
There may be some imaging error from radiation emitted at the disk boundary (being 12.5 
mm off-axis) to the detector despite the 10°collimation angle (see Fig. [I]). This effect is much 
weaker for the small resonator, where theory and experiment agreed excellently. Finally, the 
results for the large resonator with hole at ah = 47° are shown in Fig. [5] The first main 
lobe at around 120° predicted in the theoretical analysis is missing in the experimental data, 
but experiment and theory agree well for the second lobe around 200°. From the Poynting 
vector analysis, it can be concluded that the peak around 120° originates from the coupling 
region of the resonator. 

Qualitatively, all emission patterns can be interpreted in the following way. On the one 
hand, there is emission in the direction of the tangent to the resonator at the hole position. 
This best visible for the small resonator around 270° in Fig. |4| but also around 120° in for the 
large resonator with hole at ah = 47°, see Fig. |5]a-c). For the large resonator with ah = 315° 
in Fig. [5]d-f ), this contribution is masked by the waveguide and the coupling region. A naive 
picture for this tangential emission is that a fraction of the light is deflected by the scatterer 
such that it acquires a larger incidence angle on the resonator surface compared to the 
unperturbed resonator mode, and therefore is able to scatter out of the resonator. Next 
to the main lobe in tangential direction, also side lobes with higher scattering angle are 
visible. These appear as side lobes at 210° and 180° in Fig. [4] for the small resonator. 
Note that the structure at 180° almost corresponds to a radial emission from the resonator, 
i.e., almost perpendicular to the energy flow in the resonator mode. The structure around 
210° corresponds to about 45° deflection. From the simulation data, it also appears that 
similar higher deflection angles also occur towards the resonator center, even though they 
are masked in the emission pattern since they would have to pass both the resonator and 
the waveguide before they could reach the detector. 

Analogous contributions due to larger scattering angles are also visible for the large 
resonator. For example, for ah = 315° in Fig. |5j the structure at 110° corresponds to ap- 
proximately 45° deflection, and the one at 150° to approximately 90° deflection. Finally, 
the large resonator with ah = 47° also exhibits a contribution with deflection angle around 
90°, which is observed around 200° emission direction. Interestingly, in this case the corre- 
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Figure 5. Large Resonator, the hole with d = 12.5 mm sits at 315° db 5° (top row) and 47° ± 5° 
(bottom row); a) snapshot of the FDTD calculation, b) Result of the Poynting vector analysis, as 
in Fig. |4j The distance between the resonator edge and the circle on which the Poynting vectors are 
calculated is p = 4 mm. c) Experimental data (black dotted) as well as far-field radiation pattern 
predicted from the theoretical Poynting-analysis on a linear scale (red). While the second main 
lobe observed experimentally at 130° — 180° is more narrow in the theoretical prediction, good 
overall agreement is achieved, d) snapshot of the FDTD calculation, the emission direction can be 
anticipated visually, e) Result of the Poynting vector analysis, as in Fig. [4j The distance between 
the resonator edge and the circle on which the Poynting vectors are calculated is p = 4 mm. f) 
Experimental data (black dotted) as well as far-field radiation pattern predicted from the theoretical 
Poynting-analysis on a linear scale (red). While the first main lobe predicted from the theoretical 
calculation around 120° is missing in the experimental data, experiment and theory agree well for 
the second lobe at 200°, most likely due to the slightly curved waveguide in the experiment. 
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sponding emission in 45° deflection direction is not observed. One reason for this could be 
interference with light scattered from the coupling region between resonator and waveguide. 

The emission patterns for the two hole positions shown in Fig. [5] with respect to the 
coupling position substantially differ. This is in agreement with our findings in Fig. [2] The 
waveguide already perturbs the modal structure. This effect interferes with the effect of the 
hole. Both perturbations have to be taken into account in order to determine the emission 
pattern. 



V. CONCLUSION 

In conclusion, we demonstrated both experimentally by measurements in the THz domain 
and theoretically by a finite difference time domain simulation, that a finite scatterer close to 
the circumference of the whispering gallery mode resonator allows for directional outcoupling 
without ruining the quality factor. The main emission direction was tangential, originating 
at the scatterer. However, there exist side lobes with larger emission angles. Influence of 
other perturbations such as a coupling waveguide have to be taken into account in order to 
numerically determine the emission direction. With a finite difference time domain method, 
we achieved good agreement between the experimental data and the theoretical modeling in 
all investigated cases. 
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